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Abstract 


In  the  introduction  an  overview  is  given  of  the  status  of  1/f  noise  at 
the  end  of  this  five-year  contract  period.  We  indicate  that  the  1/f  noise  as 
observed  in  our  laboratory,  as  well  as  in  other  places,  quite  often  yields 
Hooge  parameters  of  order  10”^  to  10“®,  i.e.,  two  to  five  orders  less  than  a 
decade  ago.  Much  of  this  noise  can  be  seen  as  quantum  1/f  noise  which  is  the 
limiting  noise  that  can  be  observed.  Our  measurements  on  submicron  gallium 
arsenide  devices,  microwave  narrow  base  transistors,  and  gold  films  at  below 
the  Debye  temperature,  and  in  radioactive  decay  and  partition  1/f  noise,  all 
point  in  this  direction.  Over  and  above  this  limiting  noise,  additional 
1/f-like  noise  is  often  observed.  For  this,  the  standard  physical  mechanisms, 
involving  activation  energy  processes  or  tunneling  processes,  usually  apply. 
More  research  on  this  larger  1/f  noise  is  still  needed. 

In  the  section  on  experimental  work,  we  discuss  the  status  of  1/f  noise 
in  transistors  and  in  gold  films.  Also,  we  discuss  the  high  frequency  inter¬ 
valley  scattering  noise  in  gallium  arsenide  devices.  Good  agreement  with 
Monte  Carlo  simulations  is  obtained.  New  results  are  also  presented  for  noise 
in  radioactive  decay.  Both  1/f  noise  and  Lorentzian  flicker  noise  are 
observed.  The  flicker  floor  is  lower  for  lower  a-particle  energies,  in 
agreement  with  the  quantum  theory  of  1/f  noise 

In  the  theory  section  the  theory  of  quantum  1/f  noise  is  applied  to 
electron-phonon  scattering.  Explicit  results  for  the  resulting  mobility- 
fluctuation  noise  and  for  the  Hooge  parameters  involved  are  obtained.  Numeri¬ 


cal  computations  are  in  progress. 
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1.  Base  1/f  noise  versus  base  current. 


2.  Collector  1/f  noise  versus  collector  current. 

3.  Noise  of  a  "quiet"  1  pm  gold  film. 

4.  Hooge  parameter  deduced  from  Fig.  3. 

5.  Noise  of  a  0.5  pm  gold  film. 

6.  Noise  cf  a  "noisy"  1  pm  gold  film. 

7.  Diagram  o.  pulsed  high-frequency  noise  measurements. 

8.  Equivalent  noise  network. 

9.  Complete  network  including  bias. 

10.  Equivalent  noise  network. 
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12.  Diffusion  coefficient  for  GaAs  at  T  -  300K  as  a  function  of  electric 
field  strength. 

13.  Drift  velocity  of  GaAs  at  300K  as  function  of  electric  field. 

14.  Mean  free  path  in  GaAs  in  Monte  Carlo  simulation  as  function  of  electric 

field. 
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.  Introduction  and  Summary  of  Work  under  Contract 

A  long  introduction  and  sunmary  of  the  status  of  1/ f  noise  in  solids  and 
solid  state  devices  was  given  in  the  proposal  for  the  first  year  of  this 
grant,  submitted  in  1978.  At  that  time  we  indicated  that  we  were  hopeful  that 
a  breakthrough  in  the  1/f  noise  enigma  could  be  reached  if  a  number  of  proper 
experiments  were  carried  out,  together  with  new  theoretical  developments.  We 
believe  that  after  five  years  of  work  on  this  project  such  a  state  has  been 
reached.  Noise  measurements  in  transistors,  in  metal  films,  in  gallium 
arsenide  subraicron  devices,  in  radioactive  decay,  and  of  partition  noise  in 
pentodes  have  provided  much  insight  into  the  problem. 

The  basic  behavior  of  1/f  noise  is  contained  in  Hooge's  law  [1]  , 

S^fj/I2  =  <xH/fN  (l) 

where  I  is  the  average  current,  f  the  frequency,  N  the  number  of  carriers 
involved,  and  ctj  is  Hooge's  parameter.  Originally,  Hooge's  parameter  was 

-3 

thought  to  be  a  "universal"  constant  of  magnitude  2  x  10  .  Subsequent 

measurements  (1970-1984)  in  our  laboratory,  as  well  as  elsewhere  [2],  have 

Indicated,  however,  that  is  often  very  much  smaller.  Measurements  by 

Kilmer  on  microwave  transistors  indicated  values  of  order  10-7  [3].  In  metal 

films,  when  the  appropriate  correction  due  to  degeneracy  is  made,  values  of 

10-6  to  10-7  were  found  at  temperatures  below  the  Debye  temperature  [4]  .  In 

_8 

GaAs  subraicron  devices,  values  as  low  as  6  *  10  were  obtained  [5].  All  of 
this  indicates  that  a  considerable  reduction  in  1/f  noise  has  been  obtained  in 
the  last  decade.  Most  significantly,  this  low  1/f  noise  now  indicates  that 
the  fundamental  limit,  set  by  the  quantum  theory  of  1/f  noise  proposed  by 
Handel  in  1975  [6,7],  has  now  been  reached.  Thus,  much  of  present-day  1/f 


noise  work  can  be  explained  in  terms  of  quantum  1/f  noise.  Explicit 
quantitative  results  of  this  theory  for  mobility  fluctuations,  caused  by 
electron-phonon  scattering,  are,  however,  only  recently  forthcoming.  Van  der 
Ziel  and  Handel  have  recently  proposed  a  simplified  theory,  involving  both 
normal  and  Umklapp  processes  [8]  ,  which  leads  to  quantitative  results.  A 
detailed  theory,  capable  of  explaining  Hooge  parameters  of  order  10-^  or  less, 
has  been  started  under  this  grant  by  Kousik  and  Van  Vliet.  Some  preliminary 
results  appear  in  this  report  (Section  C.-I).  These  results  will  probably  be 
completed  early  in  1985.  We  will  then  be  able  to  make  a  quantitative  compari¬ 
son  between  experiments  involving  conductivity  fluctuations  and  quantum  1/f 
noise . 

Another  avenue  to  test  the  validity  of  this  theory  was  provided  by 
measuring  emission  fluctuations  of  two  kinds.  First,  as  a  nonelectric 
phenomenon,  we  proposed  in  1980  to  measure  the  1/f  fluctuations  in  alpha- 
particle  radioactive  decay.  The  results  by  Jeng  Gong  [9]  indicated  the  occur¬ 
rence  of  1/f  noise,  as  manifested  by  a  flicker  floor  in  the  Allan  variance  of 
radioactive  counting.  The  magnitude  which  was  observed  was  in  fair  agreement 
with  the  results  predicted  by  Handel's  theory.  This  work  has  also  aroused 
interest  in  the  atomic  energy  community,  such  as  in  the  Center  of  Chalk  River 
[10],  since  it  now  is  apparent  that  the  predicted  flicker  flow  sets  a  limit  in 
various  nuclear  counting  experiments. 

Secondly,  as  another  verification  of  quantum  1/f  noise,  van  der  Ziel  sug¬ 
gested  the  measurement  of  partition  1/f  noise  in  pentodes.  These  measurements 
were  performed  by  C.J.  Hsieh  [11].  A  very  good  agreement  with  theory  was 
obtained.  Altogether,  the  above  three  areas — mobility  fluctuations  (or 
diffusivity  fluctuations)  in  films  and  devices,  the  flicker  floor  in  o- 
particle  radioactive  decay,  and  partition  1/f  noise — indicate  that  quantum  1/f 
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noise  has  been  proven  to  exist,  and  most  likely  is  the  cause  of  limiting  1/f 


noise  in  various  1/f  noise  phenomena.  We  also  mention  that  such  noise  can 
explain  1/f  noise  in  phonon-assisted  hopping  processes  in  amorphous  silicon, 
as  shown  in  a  very  recent  paper  by  d'Amico  and  Van  Vliet  [12]. 

The  question  now  remains  as  to  what  causes  the  other  observed  1/f  noise 
with  rather  large  Hooge  parameter  (10-^  or  higher)  as  observed  in  certain 
devices  (long  p-n  diodes  for  example)  and  in  metal  or  semiconductor  films  at 
room  temperature  or  higher.  First  of  all,  we  suggest  that  new  experiments  be 
performed  to  establish  whether  such  noise  is  due  to  mobility  or  number 
fluctuations.  Experiments  on  transistors,  performed  both  in  Florida  and  in 
Minnesota,  have  indicated  that  the  current  dependence  for  such  noise  is 
usually  quadratic  rather  than  linear,  thus  excluding  mobility  fluctuations. 
More  experiments,  possibly  involving  a  repeat  of  Kleinpenning' s  experiments  in 
the  seventies  of  Hall  effect  and  thermo  power,  should  be  performed.  It  is  our 
surmise  that  this  larger  l/f  noise  in  all  cases  is  due  to  number 
fluctuations.  Mechanisms  involving  distributed  tunnel  ranges  or  distributed 
trapping  energies,  as  suggested  by  various  people  in  the  fifties  (see  van  der 
Ziel  [13]  and  McWhorther  [14]),  are  still  the  most  likely  candidates  to 
explain  such  phenomena.  The  more  precise  nature  of  such  physical  mechanisms 
has  in  some  cases  been  established  (see,  e.g.,  Dutta  and  Horn  [15]  for  noise 
involving  activated  vacancy  trapping  in  metal  films  at  higher  temperatures. 

In  conclusion,  we  believe  that  limiting  1/f  noise  in  semiconductor 
devices  is  presently  well  understood.  Over  and  above  this  noise,  other  1/ f— 
like  phenomena  can  occur,  the  exact  cause  of  which  has  been  pinpointed  in  some 
occasions,  while  for  others  more  study  is  definitely  needed. 

This  report  will  focus  on  the  various  experiments,  as  well  as  on 
theoretical  work,  which  were  performed  under  the  final  year  of  this  five-year 
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grant  period.  We  also  emphasize  that  the  study  of  subraicron  devices,  as 
carried  out  under  this  grant  by  Schmidt  et  al.  [5],  has  considerably  aided  <r\ 
understanding  the  importance  of  mobility-fluctuation  1/f  noise.  Great  thanks 
are  due  to  the  Cornell  National  Submicron  Resource  and  Research  Facilities 
(Drs.  E.  Wolf  and  L.  Eastman)  for  providing  the  samples  of  this  study. 

B.  Experimental  Work 

I.  Study  of  1/f  Noise  in  Transistors  (J.  Kilmer) 

The  report  of  two  years  back  (1981/1982)  contained  Kilmer's  measurements 
on  1/f  noise  in  transistors.  It  was  shown  that  by  proper  choice  of  the  source 
resistance,  either  the  base  current  noise  or  the  collector  current  noise  could 
be  isolated.  The  Hooge  parameters  for  these  noises  were  determined  and, 
especially  for  the  base  current  noise,  it  was  found  to  be  extremely  low.  It 
was  not  recognized  at  that  time  that  a  measurement  of  the  current  dependence 
of  the  noise  could  exclusively  indicate  whether  the  noise  stems  from  mobility 
fluctuations  (alias  diffusivity  fluctuations)  or  from  other  causes.  In  the 
former  case,  the  dependence  should  be  linear.  Below  we  reproduce  a  recent 
paper  by  Kilmer,  van  der  Ziel,  and  Bosnian,  reporting  such  measurements. 


Mobility  Fluctuation  1/f  Noise  in  Silicon  p+-n-p  Transistors 

by 

J.  Kilmer,  A.  van  der  Ziel  and  G.  Bosman 
Dept,  of  Electrical  Engineering,  University  of  Florida 
Gainesville,  FL,  32611,  USA 

Summary 

A  study  of  the  current  dependence  of  base  1/f  noise  and  of  collector  1/f 
noise  in  p+-n-p  transistors  shows  that  the  former  is  most  likely  of  the 
mobility  fluctuation  1/f  noise  type  and  that  the  latter  is  most  probably  not 
of  that  type.  The  current  dependence  of  1/f  noise  in  transistors  is  a 
powerful  tool  in  the  Interpretation  of  the  noise. 


In  recent  papers^’2  we  showed  that  the  low  frequency  1/f  noise  in  low- 
noise  p+-n-p  transistors  (GE  82-185)  could  be  represented  as  mobility 
fluctuation  1/f  noise.  We  found,  respectively,  for  the  Hooge  parameter  of  the 
collector  noise  =»  1.92  *  10“®  and  for  the  Hooge  parameter  of  the  base 
noise  ^  =  9.3  x  io“®,  and  attributed  these  small  values  of  ct^  to  the  decrease 
in  (Xjj  at  high  doping 


0=0:  ( — B — ) 


(1) 


where  is  the  mobility  due  to  phonon  scattering. 

But  the  theory  of  Kleinpenning  and  van  der  Ziel4  also  predicts  a 
specific  current  dependence.  For  the  collector  noise 


Sifc<f)  " 


(f)  =  2ql 


Ep 


Ep  Aft 


dp 


n  rP(0)  1 


(2) 


where  a  is  the  Hooge  parameter  tor  holes,  Ip  the  emitter  current, 


p  —  — <=>-  r  -  Ep 

Tdp  =  wB2y,2Dp  ls  the  diffuslon  cime  for  holes  through  the  base  region,  wfi  is 

the  base  width,  Dp  the  diffusion  constant  for  holes  and  P(0)  and  P(w&)  are  the 
hole  concentrations  for  unit  length  at  the  emitter  side  and  at  the  collector 
side  of  the  base,  respectively;  here  Ig.p 

Moreover,  for  electron  injection  from  base  to  emitter  we  have 


I£  is  the  collector  current. 
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(3) 


where  an  is  the  Hooge  parameter  for  electrons,  x<jn  =  wE^/2Dn,  WE  t*'e  w*-dth 
the  emitter  region,  Dn  the  diffusion  constant  for  electrons,  whereas  N(0)  and 
N(wE)  are  the  electron  concentrations  for  unit  length  at  the  base  side  of  the 
emitter  and  at  the  emitter  contact,  respectively;  here  Ign  »  Ig  is  the  base 
current.  Neglecting  the  slow  dependence  of  P(0)/P(Wg)  and  of  N(0)/N(w£)  upon 
the  current,  (f)  and  (f)  are  proportional  to  Igp  and  IEn  respectively. 

Measurements  by  Zhu  on  NEC  57807  n+-p-n  microwave  transistors^  gave  a 
different  current  dependence,  indicating  that  the  1/f  noise  in  these  devices 
was  probably  not  of  the  mobility  fluctuation  type.  It  therefore  seemed 
worthwhile  to  investigate  the  current  dependence  of  the  1/f  noise  in  GE  82-185 
transistors. 


The  results  at  l  Hz  are  shown  in  Figs.  1  and  2  for  Stfb(f)  and  Slfc(f), 
respectively.  We  see  that  Sj^b(f)  varies  as  Ig  to  the  power  (1.14  ±  0.07),  in 
agreement  with  (3),  whereas  Sj^c  varies  as  Ic  to  the  power  (1.7  ±  0.3).  We 
therefore  conclude  that  Sifb(f)  is  most  likely  caused  by  mobility-fluctuation 
1/f  noise  and  that  S^^c(f)  is  of  a  different  origin.  This  seems  to  be 
corroborated  by  the  fact  that  s^fb(f)  differs  little  from  transistor  to 
transistor,  as  we  observed  by  measuring  the  noise  of  three  different 
devices.  This  is  expected  for  fundamental  1/f  noise  mechanism  described  by 
(3).  Sifc(f)  differed  strongly  from  unit  to  unit,  as  expected  for  a  less 
fundamental  noise  mechanism  than  the  one  described  by  (2). 

The  authors  are  indebted  to  Mr.  X.C.  Zhu,  University  of  Minnesota,  for 
making  his  data  available  to  us.  The  work  was  supported  by  the  Air  Force 
Office  of  Scientific  Research  under  grant  No.  82-0226. 
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1/f  Noise  in  Gold  Metal  Films  (J.  Kilmer) 


The  previous  report  (1982-1983)  contained  Kilmer's  measurements  on  gold 
metal  films  presented  at  the  Montpellier  conference.  The  results  have  since 
been  published  in  the  open  literature  [4].  At  that  time,  however,  no  correc¬ 
tion  was  made  for  the  temperature  change  due  to  heating  by  the  d.c.  current. 
The  corrected  results  for  a  "quiet”  1  pm  device  are  shown  in  Fig.  3.  The  top 
part  of  the  figure  indicates  the  relative  noise,  S^/I2.  The  bottom  part 
gives  the  slope  y  of  the  noise  (ST/I2  =  a  J f^N  ,).  We  note  that  due  to 
heating  it  was  not  possible  to  measure  below  50K.  However,  the  trend  of  the 
temperature  dependence  of  the  noise  is  the  same  as  previously  reported.  Below 
the  Debye  temperature  (165K  for  gold)  the  noise  changes  in  character  and  is 
most  likely  caused  by  quantum  1/f  noise.  The  maximum  occurs  close  to  the 
Umklapp  temperature  (about  60K  for  gold).  At  low  temperatures  the  noise 
decreases  since  Umklapp  processes  can  no  longer  be  excited.  The  deduced  Hooge 
parameter  is  shown  in  Fig.  4.  We  note  that  the  Hooge  parameter  is  very  low. 
It  is  in  reasonable  harmony  with  the  value  computed  from  thie  theory  of  quantum 
1/ f  noise  [4]  . 

Other  devices  showed  that  the  noise  went  up  again  at  low  temperatures, 
see  Figs.  5  and  6.  These  devices  were  more  noisy  and  have  probably  not 
reached  the  quantum  1/f  noise  limit. 

A  complete  quantitative  comparison  with  theory  has  to  wait  until  G. 
Kousik's  computations  are  completed,  see  Section  C.-I. 


III.  Intervalley  Noise  in  GaAs  Mesa  Structures  (J.  Andrian) 


High  frequency  noise  at  11  MHz  and  17  MHz  was  measured  in  GaAs  diodes  of 
1.1  pra  contact  spacing  under  pulsed  conditions.  In  this  frequency  range  the 


noise  is  flat,  being  due  to  velocity-fluctuation  noise,  or  diffusion  noise. 


The  latter  is  a  function  of  the  electric  field  E,  and  is  caused  by  intervalley 
scattering.  The  noise  source  for  such  a  process  is  given  by,  see  Van  Vliet 
[16], 


SH(x,x’ ,  jw)  =*  4q2<n(x)>D(E)6(x  -  x' )  (2) 

where  D(E)  is  the  diffusion  coefficient  defined  by  the  Fourier-Laplace 
fluctuation  correlation  function 

CD 

D(E)  =  /  dte_jWt<Av(t)Av/0)>  (3) 

0 

(we  neglect  the  quantum  correction  factor,  i.e.,  hw«kT,  see  [16]).  Though  at 
high  fields  far  from  thermal  equilibrium  no  Einstein  relation  is  valid 
(contrary  to  a  claim  by  Nougier  and  Rolland  [17],  and  earlier  claims  by 
Thornber  [18]),  we  can  define  a  noise  temperature  by 


qD(E)  =  p*kTn 


(4) 


where  p*  is  the  differential  mobility.  From  integration  of  the  Langevin 
equation 


AJn  =  qnA[pn(E)E]  +  qA|>n(E)  ~]  +  H(x,t) 


it  can  now  be  shown  [19]  that  the  terminal  noise  becomes 


L  L 

Sy  -  Re  /  /  Sg(x,x' ,ju)dxdx'  “  4kT^r^Af 


(5) 


(6) 
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where  is  Che  differencial  resisCance.  Accordingly,  from  a  measurement  of 
Sy^,  Tn,  and  hence  D(E),  can  be  deCerained.  The  decalls  are  described  in  a 
paper  by  Andrian,  Bosnian,  and  Van  VlieC,  produced  under  Che  granc.  This 
paper,  which  will  shorcly  be  submicced  Co  Solid  Scace  ElecCronics  in  a 


slighdy  alCered  fora  (accounCing  for  Che  inhomogeneicy  of  E) ,  follows  below. 
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Hot  Electron  Diffusion  Noise  Associated  with  Intervalley 
Scattering  in  Very  Short  GaAs  Devices 

by 

Jean  Andrian,  Gijs  Bosraan,  Aldert  van  der  Ziel, 
and  Carolyn  Van  Vliet 

Department  of  Electrical  Engineering 
University  of  Florida,  Gainesville,  FL,  32611 


Abstract 

Noise  measurements  at  11  MHz  and  17  MHz  are  reported  for  N+N_N+  GaAs 
devices  with  an  N_  region  thickness  of  1.1  pm.  The  devices  were  fabricated  by 
molecular  beam  epitaxy  at  Cornell  University.  The  longitudinal  diffusion 
coefficient  calculated  from  these  noise  measurements  is  compared  with  that 
obtained  by  Monte  Carlo  simulation.  A  new  method  for  high-frequency  noise 
measurements  of  low  impedance  (<  IQ)  devices  is  thoroughly  discussed. 

I .  Introduction 

The  progress  in  semiconductor  device  technology  has  enabled  the  industry 
to  fabricate  devices  of  ever-decreasing  size.  In  small  devices  it  is  very 
important  to  understand  how  the  transport  coefficients,  like  mobility  and 
diffusion  coefficients,  depend  on  the  electric  field. 

Relatively  low-bias  voltages  can  easily  produce  high  electric  fields 
inside  these  small  devices.  Those  coefficients  are  needed  to  evaluate  the 
operating  conditions  of  a  device.  In  addition,  the  diffusion  coefficient 


provides  a  measure  of  Che  velocity  fluctuations  and  their  contribution  to  the 
electrical  noise  in  the  device  itself. 

In  the  case  of  GaAs,  many  authors  [  1—3  ]  have  published  theoretical 
results  (analytical  or  computer  simulations)  calculated  on  the  electric  field 
dependence  of  the  mobility  and  the  diffusion  coefficients,  but  except  for  the 
classic  work  of  Ruch  and  Kino  [4],  little  has  been  done  on  the  experimental 
study  of  high  field  transport  coefficients  [5,6]. 

In  this  paper  we  present  experimentl  results  for  the  longitudinal 
diffusion  coefficient  using  noise  measurements. 

This  method  and  the  results  will  be  described  in  Section  II.  In  Section 
III  results  obtained  using  Monte  Carlo  simulations  are  given.  The  last 
section  is  devoted  to  discussions  and  comparisons  between  our  experimental 
data  and  theoreticl  computations. 

II .  Device  Description  an  Experimental  Setup 

The  measurements  have  been  performed  on  N+N_N+  GaAs  mesa  structures.  The 

active  region  (N-  layer)  has  a  thickness  of  1.1  pm,  a  cross-sectional  area  of 

—5  o  IS  + 

7.8  *  10  vcm  and  a  Si  doping  of  10  cm  .  The  N  layers  are  needed  to 

provide  good  ohmic  contacts  at  both  ends  of  the  device.  The  block  diagram  of 

the  experimental  setup  is  shown  in  Fig.  7.  A  pulsed  bias  method  was  used  to 

avoid  excessive  Joule  heating  and  possible  destruction  of  the  device  at  high 

electric  fields.  The  pulse  generator  is  synchronized  with  a  switch  at  the 

input  of  the  power  detector  so  that  one  only  measures  the  noise  when  the  bias 

is  applied. 

Bareikis  et  al.  [6]  have  used  this  method  before;  however,  a  detailed 
description  of  the  step-up  transformer  at  the  input  of  the  amplifier  is  given 
in  this  paper.  This  network  is  very  important  because  the  device  has  a  very 


low  Impedance  (<  IQ)  compared  Co  Che  50Q  inpuC  Impedance  of  Che  amplifier. 
As  a  consequence,  Che  device  noise  vrould  drown  in  Che  amplifier  noise  wichouC 
a  sCep-up  neCwork. 

II. a.  SCep-up  cransformer 

Since  we  have  an  amplifier  wich  a  noise  equivalenC  resisCance  of  60Q 
and  a  device  wich  a  resisCance  less  chan  IQ  ,  we  need  a  tuned  circuit  to 
increase  the  resistance  seen  by  the  amplifier.  The  noise  equivalent  circuit 
is  shown  in  Fig.  8,  where  r  is  the  series  resistance  of  the  inductor,  r^  the 
device  small  signal  resistance,  and  AkTrAf  and  AkTr^Af  (k  is  the 
Boltzmann  constant)  are  the  thermal  noise  sources  associated  with  r  and  r^ 
respectively. 

The  impedance  at  the  terminals  is  given  by 


jwL  +  r  +  r 


1  -  w2CL  +  jtjdC(r+r  ) 


(1) 


If  we  operate  at  the  tuning  frequency  w  =*  - — 


/LC 


and  if  we  choose  so  that 


r  +  rd  «  |  ,  then 


Z  *  4 


C  (r+rd)  • 


(2) 


The  spectral  density  S  of  the  A.C.  open-circuited  voltage  fluctuations  at  the 
output  of  the  network  is 


S  -  4k(T0r  +  Tdrd) 


_ i_ 

jU)0C 


-—t  +  r  +  r .  +  ju)_L 

JUqC  d  0 


2 


(3) 


where  Tq  is  the  inductor  temperature  and  T^  the  device  temperature.  Carrying 


out  the  algebra  and  again  using  the  fact  that  u^CL  *  equation  (3) 

becomes : 


S  -  4k(T0r  *  V  )  i  -J—  .  (4) 

(r+rd)2 

The  condition  for  "up"  transformation  of  the  noise  equivalent  resistance  of 
the  device  is  easily  derived  from  equation  (4)  if  we  neglect  r.  Then  S 
becomes 


S  =  4kT 


d  Cr .  * 
a 


'Up"  transformation  takes  place  if 


L  .  L  .  2 

crA  >  r<i  or  c  >  cl  ■ 

a 


For  the  series  resonance  step-up  transformer ,  it  holds  that 


rd2q2 


where  Q  is  the  quality  factor  of  the  transformer.  Hence,  condition  (6)  is 
fulfilled  as  long  as  Q  >  1.  However,  a  high  Q  circuit  has  a  narrow  passband; 
therefore,  we  had  to  build  transformers  for  each  frequency  of  interest 
(11  and  17  MHz). 

II. b.  Determination  of  the  series  resistance  of  the  inductor 


If  we  replace  the  device  by  a  short  circuit,  we  obtain 


Next,  a  known  resistance  (r q)  is  used.  Equation  (4)  gives 


S2  =  4kT0  C(r+rQ) 


Equations  (8)  and  (9)  yield 


!k. !  +!o 

S2  1  +  r  • 


From  equation  (10)  one  obtains  the  value  of  the  series  resistance  r  of  the 
inductor  at  the  tuning  frequency. 

II .c.  Determination  of  the  device  resistance  with  noise  measurement 

The  actual  device  is  then  placed  in  series  with  the  inductor.  In  thermal 
equilibrium  with  both  the  device  and  the  inductor  at  room  temperature  Tn,  we 


S3  ■  4kT0  I  • 


Equations  (9)  and  (11)  give: 


r  +  fd 
r  +  rr 


Every  quantity  in  (12)  is  known  except  r^ .  From  the  noise  measurements  we 


found 


r.  *  (.65  ±  .06)Q 
a 


A  measurement  using  a  DC  ohm-meter  gave 


r,  *  .58Q  .  The  difference  between  these  two  results  can  be  attributed  to 
a 

the  skin  effect  in  the  thin  copper  wires  of  the  inductor,  which  increses  the 


resistance  of  the  wires  at  higher  frequencies. 


II. d.  Noise  measurement  with  bias  applied 


All  previous  calculations  are  valid  under  thermal  equilibrium 
conditions.  However,  we  are  interested  in  the  high  field  behavior  of  the 
noise  so  we  need  to  include  a  biasing  circuit  in  the  previous  model.  The 
circuit  is  shown  in  Fig.  9  and  the  noise  equivalent  circuit  in  Fig.  10. 


v  ■  /4kT  R  Af  ,  v,  =*  /4kT  ,r  Af  ,  and  v  =  /4kT^rAf  . 
s  Os’d  dd*  r  0 


We  take  into  account  possible  excess  noise  of  the  device  by  using  the  noise 
temperature  in  the  expression  for  v^ .  A  simple  analysis  of  the  circuit 
gives  for  the  voltage  fluctuations  at  the  A.C.  open-circuited  doutput 
terminals  at  resonance  frequency  t o 


[VvV  +  VvVi 


Since  the  noise  sources  are  not  correlated,  the  mean  cross-products  v  v  , 

s  d 

v  v  ,  and  v.v  vanish.  Hence  we  obtain 
s  r  dr 


r2(s  +S  )  R2(S  +S  ) 

d  v  V  s  V  V 

S(»„)  -  i - *  i - i— i-r 

<R8t+RSrd+t<lC)  (RS'+RSr<l+rdr) 


S  V.  V 

d  r 


2r  R  S 

.  d  s  v 

+  i ^ 


(R  r+R  r ,+r ,r) 2 
s  s  d  d 


Since  r.  <<  R_  and  r  «  R  (R  =  50Q)  ,  we  drop  the  terms  with 
u  s  “  s 

r2  and  TjR  in  the  numerator.  Equation  (14)  becomes, 
d  d  s 

<S  +S  .  ) 


$ 


with  Sy  =  4kTQf  and  Sy  =*  4kT(jrd.  Equation  (15)  can  be  written  as 


S(wn)  -  4k  £  - - -  (T.r+T.r.)  . 

0  C  rKl+f-)2  0  dd 


(16) 


Equations  (8)  and  (16)  yield 


r2(l+  — 12T 
S1  d^  rd  ° 
S"=  ^(T0r+Tdrd) 


(17) 


or 


(r+r,)% 


T  r  +  T  r  =  — 
dd  x0  Sl 


(18) 


We  define  the  noise  equivalent  resistance  of  the  device  as 


T  r 
d  d 


(19) 


From  equation  (18)  we  obtain 


s  (r+rd)2 

Rn  .  — - r 


(20) 


In  Figure  11  Rr  is  plotted  as  function  of  the  applied  voltage  for  tq  *  300K 


w, 


and  f. 


0 


11.0  MHz  and  fn  *  17.0  MHz.  We  note  that  for  bias  voltages  up 


‘0  2*  . . .  0 

to  0.2V  the  measured  value  of  Rr  is  equal  to  the  small  signal  resistance  of 

the  device  rd,  within  the  experimental  error.  For  bias  voltages  above  0.2V, 

Rn  increases  sharply  with  the  voltage  applied. 

We  define  the  noise  temperature  Tj  by  [7]  as  qD(E)  -  kT,p,  where  E  is 
— — —  a  a  a 

the  mean  d.c.  field  in  the  diode  and  is  the  differential  mobility.  It 
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can  then  be  shown  that  the  terminal  voltage  noise  is  4kT  r.Af;  then 

n  d 

T  =  T,.  We  now  find  for  the  diffusion  coefficient 
n  d 


D(E) 


R 


(20a) 


where  we  replaced  by  ~  T  (see  equation  (19).  Measurements  of  the  white 

d 

noise  and  the  current-voltage  characteristic  enable  us  to  determine  experimen¬ 
tally  the  diffusion  coefficient  as  function  of  the  applied  electric  field.  In 
Figure  12  we  plotted  the  diffusivity  D(E)  at  300K.  as  following  from  Figure 
11.  In  order  to  have  some  means  of  comparison  with  the  theory,  we  performed 
computer  simulations  which  we  describe  in  the  next  section. 

III.  Monte  Carlo  Simulation 

To  predict  the  high-field  transport  properties  of  a  given  semiconductor 
material,  one  has  to  solve  the  Boltzmann  transport  equation.  This  equation  is 
very  difficult  to  solve  analytically  unless  drastic  simplifications  are 
introduced.  Such  simplifications,  however,  may  alter  very  much  the  physical 
meaning  of  the  solution. 

To  circumvent  this  difficulty,  numerical  techniques  have  been  developed 
to  obtain  accurate  results.  One  of  these  is  the  Monte  Carlo  simulation 
technique. 

The  principle  of  this  method  is  to  simulate  on  a  computer  the  motion  of 
one  electron  in  the  k-space  through  a  large  number  of  scattering  processes. 
Excellent  papers  have  been  written  on  this  subject  by  Fawcett  et  al.  [1]  and 
more  recently  by  Boardman  [2].  For  more  details  we  refer  the  reader  to  these 
articles . 

Our  simulation  is  based  on  the  method  developed  by  Boardman  [4].  To 
explain  the  noise  behavior  of  our  device,  we  need  to  know  the  drift  velocity 


23 


and  the  diffusion  coefficient  of  the  electrons  as  a  function  of  electric  field 
strength.  The  drift  velocity  follows  from  the  average  displacement  of  an 
electron  in  a  sufficient  long-time  interval.  The  diffusion  coefficient  is 
calculated  as  follows: 

D  =  . 5(<V2>  -  <V>2)T  (21) 

where  T  is  the  sampling  time,  V  is  the  average  electron  velocity  for  a  sample 
over  this  time.  The  average  is  performed  over  an  ensemble.  In  the  computer 
simulation,  however,  we  use  time  averaging.  This  is  allowed  because  we  deal 
with  a  stationary  process,  so  the  time  average  is  equal  to  the  ensemble 
average  of  a  physical  quantity.  D(E)  and  V(E)  are  shown  in  Fig.  12  and  Fig. 
13  respectively. 

IV.  Discussion 

We  obtained  good  agreements  between  the  experimental  and  the  computer 
simulation  data.  The  intervalley  deformation  potential  is  the  only  parameter 
which  we  can  adjust  in  our  Monte  Carlo  simulation;  all  other  parameters  are 
generally  well  known.  Pozela  et  al.  [3]  used  a  value  of  10®eV/cm  to  match 
their  simulation  with  the  existing  experimental  data.  In  our  case  we  use  the 
value  of  10^eV/cra.  This  value  has  also  been  used  by  Fawcett  et  al.  [1].  We 
do  not  claim  that  this  is  an  indirect  determination  of  the  intervalley 
deformation  potential  since  more  measurements  in  a  wider  range  of  temperature 
would  be  needed  to  obtain  an  accurate  value  for  this  parameter. 

In  Figure  14  we  plotted  the  mean  free  path  of  an  electron  as  function  of 
the  applied  field,  calculated  using  the  Monte  Carlo  program  for  T  =  300K.  In 
the  range  of  electric  fields  of  our  measurements,  the  mean  free  path  is  small 
compared  to  the  length  of  the  device.  This  assures  that  the  measured  data 


A  difficulty  one  might  encounter  in  the  measurement  of  high  field 
properties  of  GaAs  is  the  possible  occurrence  of  traveling  dipole  domains. 
These  domains  would  make  the  electric  field  inside  the  device  dependent  on 
position  and  time.  However,  according  to  Kroemer  [8],  these  dipole  domains 
are  formed  only  when  the  produce  n  (doping)  by  A  (length)  exceeds 
2  *  lO^cm  .  This  is  not  the  case  in  our  device  which  has  a  doping  of 
lO^cm-^  and  a  length  of  1.1  x  10“^cm. 

Inside  the  device  the  electron  experiences  a  higher  effective  electric 
field  than  that  calculated  from  the  applied  voltage  (V/A).  This  is  due  to  the 
potential  profile  which  exhibits  a  minimum  whose  position  depends  on  the 
applied  voltage.  For  high-bias  voltages  this  potential  minimum  Is  located 
close  to  the  cathode.  The  actual  electric  field  Is  the  gradientof  this 
profile,  and  its  value  at  the  anode  side  of  the  potential  minimum  is  probabl; 
somewhat  larger  than  V/A.  As  a  consequence,  the  experimental  curve  has  to  be 
shifted  slightly  to  higher  electric  field  values. 

We  made  sure  that  we  measured  velocity  fluctuation  noise  at  11.0  and  17.0 
MHz  by  measuring  also  the  noise  at  lower  frequencies.  There  we  observed  a 
small  amount  of  1/f  noise,  characterized  by  a  Sooge  parameter  of  4  *  10-^.  No 
other  excess  noise  sources  were  detected. 
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IV.  Noise  in  Radioactive  Decay  (G.  Kousik) 

New  measurements  on  counting  statistics  of  a-particle  radioactive  decay 
in  gjAm,  g^Pu,  95cm,  and  ^Gd  were  made,  in  order  to  see  whether  deviations 
from  Poisson  statistics  occurred.  We  measured  up  to  intervals  of  2000 

minutes;  however,  since  the  counting  equipment  was  only  available  to  us  for  a 
limited  period,  not  many  long-time  intervals  could  be  measured.  So  the 

accuracy  of  the  Allan  variance  for  these  intervals  is  rather  poor.  Neverthe¬ 
less,  it  was  clear  that  the  noise  is  not  purely  Poissonian.  Both  a  Lorentzlan 
and  a  1/f  flicker  noise  component  could  be  detected  in  these  radioactive 
sources.  The  theory  and  experiment  of  these  measurements  are  described  below. 

Alpha  Particle  Counting  Experiment 
G.  Kousik,  W.H.  Ellis,  G.  Bosman  and  C.M.  Van  Vllet 

If  we  assume  the  noise  spectrum  of  a-partlcle  emission  to  consist  of 
white  noise,  1/f  noise  and  a  Lorentzlan  flicker  noise,  then 

S  (w)  =  2mn  +  2uC  1/M  +  4B  - - -  .  (1) 

m  u  2.2 

or  +  ur 


Using  the  Allan  variance  theorem  (1),  we  have  for  the  relative  Allan  variance 


R(T)  =  +  —  log  2  + - - —  [4e-0tT  -  e'2aT  +  2aT  -  3]  .  (2) 

m0T  m2  m2T2a2 

Here  is  the  counting  rate,  T  is  the  time  interval  of  counting,  and  B  and 
a  are  constants.  For  T<<a  \ 


,  -aT  -2aT  _  „  ,  „  q  ■> 

4e  -  e  +  2aT  -  3  =  2/3  a3T3  . 


Therefore 

R(T)  -  “  +  —  log  2  +  —  aT  .  (3) 

0  -8  H 

According  to  Handel's  quantum  1/f  noise  theory,  the  flicker  floor  2C  log  2/m2 


is  given  by 


-2-C-  -1-0-8--  =  33.28  x  l(f7  log  2  x  -  . 
2  < 


Here  E  is  the  energy  of  the  a-particle  in  MeV.  If  the  source  has  no  Lorentzian 
flicker  noise,  then  for  large  T  the  white  noise  term  is  much  smaller  than  the 
flicker  floor  which  is  proportional  to  E. 


R(T)  -*■  2C  log  2/m2  (for  large  T)  . 


Experiment 

Here  two  sources  with  four  different  energy  a-particles  were  used.  The 
experimental  setup  is  shown  in  Fig.  15.  To  count  adjacent  intervals  a  program 
was  set  up  in  the  autosequence  mode  of  the  ND66 ,  Multichannel  Analyzer.  This 
also  contained  a  step  for  printing  out  the  total  number  of  counts  in  each  peak 
after  several  interval  counts. 


5.155 


2.41  x  10' 


5.486 

5.805 

3.183 


Since  the  average  counting  rate  mg  will  be  different  for  different  time  Inter¬ 
vals  T,  all  of  the  relative  Allan  variances  were  normalized  as  follows: 


VT>  -  R(I>  -  +  SJr 


where  mg  is  chosen  to  the  nearest  1000  or  500  counts.  Here  mg  are  as  follows: 


94Pu 

18500 

per  minute 

g  5  Am 

13000 

«•  M 

96Cm 

10000 

•  •  M 

64^ 

18000 

•  •  tt 

Figs.  16  and  17  show  R^T)  vs.  1/T  for  these  sources. 

The  Pu  and  Cm  sources  show  a  trend  towards  a  flick  r  floor  with  that  of 
Cm  being  greater  than  that  of  Pu.  This  is  in  agreement  with  Handel's  theory 
in  that  the  flicker  floor  must  be  proportional  to  the  energy  of  the 
a-particle.  The  Am  source  shows  no  flicker  floor  but  a  Lorentzian  for  the 
longer  time  intervals.  The  Gd  peak  was  not  sharp  because  this  was  a  thick 
source  and  that  led  to  an  energy  loss  for  the  a-particles,  broadening  the  peak 
on  the  lower  energy  side.  This  source  also  shows  no  distinct  flicker  floor, 
but  a  tendency  towards  a  slow  Lorentzian  flicker  noise. 
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C.  Theoretical  Work  (G.  Kousik) 


Below  we  Indicate  computations  of  quantum  1/f  noise  for  various 
scattering  mechanisms.  For  that  purpose  we  must  first  generalize  Handel's 
expression  for  noise  in  a  scattered  beam,  so  as  to  Include  directional 
calculations. 

Handel's  original  expression  for  the  1/f  noise  in  a  particle  current  K, 
which  is  scattered  in  a  direction  with  velocity  change  Av  ,  was 


Sr/K2  =  2oA/f 

where  a  is  the  fine  structure  constant  (137)“^, 


(7) 


A 


8  |  A v  |  2  _  8h2  1  Ak|  2 

3lt  c2  =  311  m*2c2 


(8) 


and  where  m*  is  the  effective  mass  and  c  the  velocity  of  light.  In  phonon 

scattering  we  need  the  correlation  between  the  particle  current  w^  , 

scattered  from  k  to  k'  ,  and  the  particle  current  w  „  ,  scattered  from 
~L  ~  -2- 
k^  to  k"  .  A  logical  Ansatz  for  this  correlation  is 


2aA 


12 


f  "k^k*  wk^k' 


(9) 


We  took  here  the  cross  correlation  to  be  bilinear  in  both  fluxes,  analogous  to 
van  der  Ziel's  expression  for  quantum  partition  1/f  noise,  see  [11].  The 
cross  coupling  constant  should  be  given  by 


k'  =  k^  ±  q'  (+  sign  phonon  absorption)  (11a) 

k”  =  k.  ±  q'  (-  sign  phonon  emission)  (lib) 


with  k  representing  the  electron  wave  vector  and  qthe  phonon  wave  vector.  The 
above  holds  for  normal  processes.  For  Umklapp  processes, 


k'  "  k  =  ±  q'  +  g'  =  g’  (12a) 

k"  -  ^  =  ±  q”  +  g"  =  g"  (12b) 


where  g'  and  g”  are  reciprocal  lattice  vectors.  Thus 


12 .Umklapp 


8ft2 

3it 


m* 


(13) 


For  the  smallest  Umklapp  process  (adjacent  Brillouin  zones)  g'  *  g"  *  2n/am in 
=  2n/ ag,  and  we  have 


A12 ,Um 


32 n  di2 

3  (m*aQ)2 


(14) 


We  now  have  complete  expressions  for  the  scattering  correlation  spectrum 
S  where  the  w' s  are  the  rates  occurring  in  the  matrix  elements  of 

U  M 

kxk*  ’  k2k” 

the  Boltzmann  equation.  Thus  an  exact  theory  for  the  fluctuations  in  mobility 
can  be  constructed.  This  has  so  far  been  carried  out  for  acoustical  and 


optical  phonon  scattering,  involving  only  normal  processes.  The  Umklapp 
processes  will  be  considered  in  the  near  future. 


Below  we  give  a  sketch  of  two  detailed  calculations,  to  be  published  at  a 


later  date.* 


I.  Quantum  1/f  Noise  Due  to  Acoustical  Phonon  Scattering 
(G.  Kousik,  C.M.  Van  Vllet  and  G.  Bosman) 

We  start  with 


n 


e  <v2t> 

m*  .  ’ 

<vz> 


where  the  averaging  is  over  a  Brillouln  zone  in  k  space.  This  yields  for  the 
mobility-fluctuation  noise, 

S  (f) 

^n 


We  have  from  the  solution  of  the  Boltzmann  equation  (see  Figure  18), 


cos  X* 
cos  X 


em 


+  wabs(i 

wkk* 


cos  X  abs 


SAt(£  )Ax((fT.  ) 

-1  -2 


*2<£,  )*2<£k)  l  l 

h  -2  k’  k" 


cos  x' ■> 
cos  %1J 


(1 


cos 

cos 


X’ 


)  +  sab,l  (1 


w 


abs 

k^k" 


cos  X'  \  _  cos  x"-) 

cos  Xj/  k  cos  X2J 


Here  the  noise  spectrum 
contributions  from  absorption 


SAT(£k  )At(£,  )  is 

and  emission,  we  write 


written 


as 


a 


sum  of 


*We  note  that  in  all  cases  Hooge's  law  S^/p2  =*  o/fH  follows  completely 
from  the  assumptions  (9)  and  (10)  of  the  'theory.  It  is  particularly 
gratifying  that  the  factor  N  appears  automatically  in  the  denominator  as  a 
consequence  of  k-space  integration.  Explicit  results  for  the  Hooge  parameter 
are  deduced. 


s7<VV-k'’-k"> 


12  ,-,abs  nabs 


«5  rk  k_  k'  k”  ^  =  _ —  Uaua  W  & 

W  'h  »-  }  f  Wk1,k’  k^k”  ^,^2 

«Tk.  -  ijr  l^q)l2[F(q)  +  l]6(£k.  -  £k  +  ««  )  ;  v  -  k-  3 
q  .  -  a 

wkbk,  “  ^  ir  l5^q>!2F<q)6<£v  ~  £k  j  If’  -  ^  +  q 
q  - 

F(q)  =  [expC-Kw  /kgT)  -  1]  1  ;  |<^'(q)|2  =  (fic2/2  pVw  )  |  q  |  2 


(i  -  ^ssjsL)  a±iCOsp 
v  cos  x  ‘  k  v 


tan  x  sl-n  P  cos 


r+  for  emission  i 
*-  for  absorption* 


Let  us  assume  that 


□ 


“*12  - - T7  ’ 

3icra*2c2 

see  the  introduction  to  this  section.  Now  consider  the  emission  terra: 

i  -  X  X  {(|^}2|^q,>l2r^q">l2lF(q,>  +  WC k,  -  £k  +^) 

k’k"  3’  ,3"  -  -1 

V.^-q,6(^"  ‘  ^  +'KV)  x  (cos  Pl  “  tan  XlSin  picos  *i)] 

tl  ft  tV/  ^ 

[t~  (cos  P,  -  tan  x2sin  P2cos  4>2 )  ]  ( - 7-7)  7  t  q’q”}  • 


Summing  over  k'  and  k”  we  get, 


2c2'  f  W 


'em  '  tl1'2  +  U[F(,-)  +  1] 
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X  (cos  p^  -  tan  x^  sin  p^  cos  ^Xcos  ~  tan  ^  sin  P2  COS  ^2^q'  q"  * 

6 

I  ,  AuS._vL-/  /  d3q’d  v  k’VJ|?-(q')|2|^rq")|2 

°»  *  3m«2c2  fklk2  (8«3)2 

[F(q’)  +  l|[F(q")  +  1)6(6  *  £y*  *"  J  *  S(  £  ^  £^+  «*  jJ 

(cos  Px  cos  P2  +  tan  xL  tan  x2  sln  ^  sln  P2  cos  ^  cos  $2 
-  cos  Px  tan  X2  sin  P2  cos  $2  -  cos  P2  tan  xL  sin  pL  cos  4^)}  . 

We  have 

//  d3q'd3q"  =  //  q' 2dq' q"2dq"  //  d  cos  P^  d  cos  P2  //  d^d^  . 

Let  us  carry  out  the  //  d^di^  .  Writing  only  the  last  parenthesis  we  have 
2u  2n 

/  /  d^.d#  [cos  p.  cos  P2  +  tan  %l  tan  sin  PL  sin  P2  cos  ^  cos  <|»2 

0  0 

-  cos  P2  tan  X2  3in  P2  cos  ♦ i  ~  cos  $2  Can  X1  Sln  ^1  C°S  ^  * 

Since 

2n  2« 

/  /  d<|>  d$  (cos  4»  cos  4>  -  cos  #  -  cos  <t> .)  «  0  , 

Art  12  1  *•  *• 


33 


we  are  left  with 


2n  2n 

/  /  d<t>  d<|>  =  4 it2  . 

0  0 

Therefore 
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/  dq'd  cos  {q^l^tq')  |2[F(q’)  +  1 1  k^-q*  “  +  ^q.  )cos  PLJ 


x  /  dq"d  cos  ^  {q"4  iS^q")  I  2[F(q")  +  1]6(£  ^_q..  -  )cos  p 


Let  us  write 
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where  we  used  w  '  ■  c.q'  . 

q  x 


Define 


-Kg'  _  *kl 
2  m*  m* 
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and 


Q'  -  ifq"/2m*  . 


With  this  change  of  variables, 


dq'dcos  p 
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d(Q",y®) 

and 


In  the  integrals  we  get  6(-Kq'y*)  and  e6('fiq“y2)  >  which  can  be  written  as 

1  G  1  G 

6(yx)  and  6(y2)  respectively.  Therefore 
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Since  we  have  a  6(ya  )  ,  integration  is  confined  to  a  strip  at  ya  =  0 

1  )  L  A.  f  L 

along  Q'  and  Q"  (see  Figure  19).  Thus  the  limits  are 
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The  absorption  terra  is 
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Summing  over  k'  and  k”  we  get 
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X 


{r3tq")|2F(q”)6(^+a..  -  -  4uy  ) ( cos  P2  +  tan  ^sin  P2cos  <(>2)}q" 


2n 

As  in  the  emission  case  only  the  //  cos  6  cos  |3  d<{>  d<(>  term  survives 

0  1  2  12 

Therefore , 
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Consider  only  the  first  integral: 
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The  numerator  has  a  form 
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II.  Quantum  I/f  Noise  Due  to  Optical  Phonon  Scattering 
(G.  Kousik,  C.M.  Van  Vliet  and  G.  Bosaan) 
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The  relaxation  time  is 
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From  Handel’s  theory 
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both  independent  of  q. 


Call  |^"( q)|2  =  Cg  and  F(q)  =  F.  Consider  the  emission  term  in  the  sum 
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/  d<J>1d4>2  =*  4n' 


The  two  integrals  in  Iem  look  alike.  Consider 
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Integrate  along  a  strip  (see  Figure  22) 
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Now  consider  the  absorption  term 
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